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We show, by means of a simple compilation of the available experimen- 
tal results, that the pj_-spectra obtained at RHIC and SPS are strikingly 
similar up to p± ~ 1.5 — 2 GeV. In fact, the discrepancies between various 
experimental groups working at the same experiment are of the same size as 
the differences between RHIC and SPS. Our observation is complementary 
to the well known fact of the equality of the measured i? s idc and i? ou t HBT 
radii at RHIC and SPS. In essence, it points out that the transverse size of 
the firecylinder and the strength of the transverse flow are not significantly 
changed between SPS and RHIC. This suggests that a saturation mecha- 
nism is effective already at SPS. We also point out that the dominance of 
protons over ir + at large p± can be seen not only in the RHIC data, but 
already in the SPS data. 

PACS numbers: 25.75.-q, 25.75.Dw, 25.75.Ld 

1. Introduction 

The pertinent question in the field of relativistic heavy-ion collisions is 
whether the physics at RHIC is qualitatively different from the physics at 
SPS. In our opinion, the available experimental results hint that this is not 
the case, at least for soft processes. 

In this paper we compile the transverse-momentum spectra of hadrons 
measured by various groups at SPS and RHIC. Surprisingly, to our knowl- 
edge such a study has not been presented before. We use the data of NA44 
0], NA49 §, PHENIX §, and STAR [|,|, §, and show that there exist 
discrepancies between the NA44 and the NA49 data, as well as between the 
PHENIX and the STAR data. In fact, these discrepancies are of the same 
size as the differences between SPS and RHIC. More precisely, within the 
experimental uncertainties, which are quite large, one finds that the slopes 
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of the pj_-spectra at RHIC are compatible with those at SPS. This observa- 
tion is complementary to the well known fact of the very weak beam-energy 
dependence of the transverse HBT radii. Similarity of both the transverse 
size and the pj_-spectra of hadrons indicates that the amount of the trans- 
verse flow cannot be significantly different at the two considered collision 
energies. 

The transverse HBT radii R s ide and R ou t measured [f7|, ||, ||] in the first 
run of RHIC are very close to those measured in heavy-ion collisions at 
smaller beam energies. Only the longitudinal radius Ri ong exhibits a mono- 
tonic growth with y / sjv r /v (for a compilation of the data at different ener- 
gies see, e.g., Ref. 0). The weak energy dependence of the transverse 
radii is surprising, since the RHIC beam energy, ^snn = 130 GeV, is al- 
most one order of magnitude larger than the SPS energy, y / ijv r 5v =17 GeV, 
and one would naively expect that much larger hadronic systems were pro- 
duced at RHIC. One would also expect a longer lifetime of the hadronic 
fireball formed at RHIC, which should be reflected in a longer emission 
times of pions. This effect is quantified by the measurement of the ra- 
tio R ou t/ Rsidei which is expected to be much larger than unity for long 
emission duration [10]. The experimental measurements indicate, however, 
that Rout /Rside is compatible with unity in the whole range of the studied 
transverse- momentum range (0.2 < kj- < 1.0 GeV). This fact is another 
puzzle delivered by the analysis of the RHIC data. 

On the other hand, the first measurements at RHIC showed that the 
pseudorapidity densities of charged particles are higher than those observed 
at SPS. Can this effect be reconciled with practically unchanged transverse 
radii? For the most central collisions, PHOBOS communicated the value 
{N ch } = 555±12±35 Q, BRAHMS the value (N ch ) = 549±1±35, whereas 
PHENIX obtained (N ch ) = 622 ± 1 ± 41 p|. Normalizing per participant 
pair yields (N ch ) /(0.5(N part )) = 3.2 for both PHOBOS and BRAHMS, and 
(N c h) I '(0.5 /N part)) = 3.6 for PHENIX. These numbers may be compared to 
the NA49 result, 1.9 @, 0, and the WA98 result, 2.6 @. We can see that 
the multiplicity increases by about 50% when we move from SPS to RHIC. 
A simple geometric scaling suggests that the transverse radius increases, 
correspondingly, as a square root, i.e., by about 20%. Thus, the observed 
increase of the multiplicity translates to moderately small increase of the 
transverse radii. Clearly, the difference between the NA49 multiplicity and 
the WA98 multiplicity, as well as the errors of each particular experiment 
lead to the uncertainty in the determination of the geometric parameters. 

The weak dependence of R s %de from SPS to RHIC means that the trans- 
verse size of the firecylinder changes very little. If this is the fact, than the 
amount of the transverse hydrodynamic flow should also be similar, since 
it is difficult to imagine that a much stronger flow would lead to the same 
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transverse size at freeze-out. To the contrary, the STAR results for the 
P±-spectra have been interpreted as an indication of a much stronger flow 
at RHIC. Below we redo the simple analysis of the hadronic spectra and 
show that the errors in the flow parameter are very large, such that one 
cannot definitely conclude that the flow is stronger at RHIC. Moreover, the 
PHENIX data suggest a much lower flow than STAR, such that it becomes 
compatible with the flow at SPS, as obtained from NA44 and NA49. 



2. Compilation of the p^-spectra measured at SPS and RHIC 

We begin by just displaying the experimental results from various groups 
in a single plot. In Fig. 1 we plot the pj_-spectra of tt~ , K~ , p, and A. 
Similarly, in Fig. 2 we show the pj_-spectra of tt + , K + , p, and A. The 
collected data come from the measurements done at midrapidity for the 
most central events ffl, 0, pi 0, S, pi . A striking feature of Fig. 1 is a very 
impressive agreement between different experiments for the pions. In the 
range 0.4 < p± < 1.0 GeV, the data from NA49 coincide with the data from 
PHENIX and STAR, whereas the data from NA44 show the same slope 
with a slightly smaller normalization due to a different centrality choice 
(see the discussion below). A very similar pj_-dependence is also seen in the 
spectra of kaons, p, and A measured by different experiments. In this case, 
however, a different normalization of the spectra between SPS and RHIC 
data is clearly seen. In Fig. 2 one can clearly see similarities in the shapes 
of the pion as well as the kaon spectra measured by NA44 and PHENIX. 



If the plots for mesons and baryons were overlaid (c/. [15, |16[), one 
could see that the PHENIX data have the property that the pj_-spectra of 
7r + and p cross around p± = 2 GeV, such that there are more protons than 
7r + at large momenta. We wish to stress that the same phenomenon can be 
seen already at SPS, where the protons dominate ir + at p± > 1 GeV. The 
lower value of the crossing point reflects a much higher proton density at 
SPS compared to RHIC. 

At RHIC, the "anomalous" (3|, [17]] feature of the spectra is the dom- 
inance of p over ir~ at pj_ > 2 GeV. One could speculate that a similar 
behavior might occur already at SPS, if the measurements had beed carried 
out to sufficiently high momenta. Indeed, when the NA44 data are ex- 
trapolated with simple exponential functions, then one finds that a crossing 
occurs around p± = 2.5 GeV. Since the exponential fits may not work over 
the large range in p±, this phenomenon remains a speculation till verified 
exp er iment ally. 

In order to see the similarities between the data even more vividly, we 
scale the normalization of the spectra in the following way: 
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In the NA49 data we undo the corrections for the feeding of protons 
and antiprotons from weak decays. According to Ref. this correc- 
tion is about 30%, thus we divide the NA49 data for p and p by the 
factor 0.7. 

The most central data from NA44 correspond to the average impact 
parameter (6) = 5 fm, whereas the most central NA49 data correspond 
to (b) = 2 fm Q. This difference explains a smaller normalization of 
the NA44 data. We correct for the centrality choice of NA44, multi- 
plying all NA44 spectra by an educated-guess factor of 1.5. 

Inspired by the success of the thermal approach, we scale the spectra 
by the factor exp[— (BfiB + Sfis + hfJ-^/T], where T is the freeze-out 
temperature, B, S, and I3 are the baryon number, strangeness, and 
the third component of the isospin of the particle, respectively, and 
Hb, /As, and fii are the corresponding chemical potentials. At SPS 
(Pb + Pb, = 17 GeV) m 



T = 164 MeV, \i B = 229 MeV, fi s = 54 MeV, m = -7 MeV, 

(1) 

while at RHIC (Au + Au, ^vlv = 130 GeV) |l|] 

T = 165 MeV, y, B = 41 MeV, y, s = 9 MeV, m = -1 MeV. (2) 

Since in the thermal model the (original) spectra are proportional 
to the factor exp[(£?/i£ + Sfis + IsHi)/T] (for the Boltzmann statis- 
tics, which works very well), our rescaling approximately removes the 
effects of different chemical potentials at RHIC and SPS. Due to reso- 
nance decays, even within the thermal model the scaling is not exact, 
but approximate, f 

The described rescaling is not necessary for our analysis, and has no 
impact on the conclusions. However, it is useful, since it brings the spectra 
closer and makes the eye-ball comparison easier. The results are shown in 
Figs. 3 and 4. For the pions, the four measurements agree very well in the 
overlap region. For the kaons, the NA49 data have a visibly larger slope, 
while the three other sets of data overlap. The measurements for protons 
and p agree as well. The spectra of A and A have smaller slopes in the 
STAR measurement compared to NA49. One should bare in mind that 
the statistical errors are typically of the order of a few percent, and the 
systematic errors are around 10%. Hence, within the experimental errors, 



1 For instance, a larger baryon chemical potential leads to a larger pion yield from the 
process A — > nN, hence the secondary pions are sensitive to the baryon chemical 
potential. 
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there are no significant differences between SPS and RHIC which could 
indicate different physics. Moreover, the small difference between NA44 
and NA49, and between PHENIX and STAR, which one may observe in 
Figs. 3 and 4 is of the same magnitude as the difference between SPS and 
RHIC. 



3. Exponential fits 

In this Section we wish to quantify the eye-ball observations from Figs. 3 
and 4. The most popular method is to fit the spectra to the exponential 
function, and compare the parameters. In fact, this is the most economic 
and common way of presenting the data, and the differences in the slope 
parameters of particles of different masses are interpreted as the signature 
of the transverse flow 21]. Hence, we fit the functionP] 



d 2 N 



m±dm±dy 



= Aexpi-m^/X) (3) 

y=0 



to each rescaled spectrum, independently for different experimental groups^] 
The method is as follows: we constrain the fit to p± < 2GeV (this has 
relevance for the PHENIX data, and for A and A measured by STAR). For 
simplicity, we assume a 15% error on every point. This is in the ballpark 
of the errors given by the experimental groups. Then, the x 2 function is 
minimized with respect to the A and A parameters. 

The results are shown in Fig. 5. The first feature to notice are large cor- 
relations between A and A, which lead to sizeable errors in these parameters. 
The optimum values for A, denoted in Fig. 5 by crosses, and the errors are 
listed in Table |l]. While the optimum values agree with the estimates made 
in other papers [|], [| ||, f|], our errors are significantly larger for the case of 
STAR and NA49. The physically relevant parameter is the inverse slope, 
A, while the norm parameter, A, carries the ambiguities described above 
(rescaling). Hence, if the results of the two data sets are to be consistent, 
then the ranges of A (and not necessarily A) should overlap. We note from 
Fig. 5 that this is basically the case for NA44 and NA49, except perhaps 
the case of tt + . The STAR and PHENIX data are much less consistent, 
especially for the case of tt~ . For kaons, protons, and p, the PHENIX data 
agree within errors with NA44 and NA49. Some caution is needed in the 



2 Other forms used in the literature, differing by the power of m±, lead to similar 
qualitative conclusions. 

3 We note that the slope parameters, A, fitted in many papers, depend strongly on 
the choice of the range in p± (see, e.g., the discussion in P2]). Here we use the 
available data ranges with p± < 2 GeV. Because of the lack of strict thermodynamic 
interpretation, we refrain from calling A the temperature of the spectrum. 
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Table 1. Our values for the inverse-slope parameters, A, in units of MeV, fitted 
to mid-rapidity data from the most-central collisions of Pb + Pb at ^/snn = 17 
GeV (NA44 §, NA49 §) and Au + Au at = 130 GeV (PHENIX §, STAR 

[§ |, U). The fit includes data points up to 2 GeV. The NA49 data, and the STAR 
data for 7r~ and K~~ are preliminary. 



interpretation of Fig. 5, since the fits are made over different ranges in p±, 
and the form used for the fit has no sound physical ground. Yet, there are 
indications for discrepancies between STAR and PHENIX. Also, it is clear 
that only the STAR data suggest larger values of A (albeit with large er- 
rors), which would mean larger transverse flow. If PHENIX data are used, 
no such conclusion comes out, since the PHENIX data are basically consis- 
tent with SPS. Note also that the quality of the fit, indicated by the value 
of x 2 a t the optimum, is worst for the case of pions from PHENIX. This 
reflects the inapplicability of the simple exponential parameterization over 
the large range in p±. 

The final remark of this Section concerns the multiplicities of particles, 



dN 
dy 



y=o 



/ dp±p± 
Jo 



d 2 N 



p±dp ± dy 



(4) 

y=0 



The integrand in Eq. (|j) peaks at p± ~ 0.25 GeV for the pions, and 
0.6 GeV for the protons. Therefore much of the strength comes from the 
p± region not covered by the data. The necessary extrapolation brings 
systematic uncertainties to dN/dy and may be a source of discrepancies 
between various quoted numbers. 

4. A model calculation 

The estimates of the transverse flow must be based on model calcula- 
tions which include this effect, as well as other potentially important phys- 
ical effects. In this Section we apply the single- freeze- out model, which 
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has been used successfully to describe the pj_-spectra at RHIC. The model 



has been described in detail in of Refs. [15, 16, 22], and here we only list 
its basic features: i) simultaneous chemical and thermal freeze-out of the 
hadronic matter, ii) inclusion of all hadronic resonances, and hi) a sim- 
ple parametrization of the freeze-out hypersurface, which is defined by the 
condition r = \Jt 2 — x 2 — y 2 — z 2 =const. The hydrodynamic flow on the 
freeze-out hypersurface is taken in the form resembling the Hubble law, i.e., 
u n _ X n J T- xhat way both the longitudinal and transverse flows are built 
in. The single-freeze-out model has two thermodynamic parameters (tem- 
perature and the baryon chemical potential) which are fixed by the global 
fit to the relative particle yields (c/. Eqs. (]lj,@)). The extra two parameters 
(r and the transverse radius of the firecylinder /0 max ) determine the over- 
all normalization and the shape of the spectra. In our present calculation 
the thermodynamic parameters are the same as those used previously in 
Sect. 2. Q Then, the expansion parameters r and p max are fixed separately 
to the data of each experimental group. They determine, in each particular 
case, the maximal transverse flow, given by the model formula 

amax _ Pmax /_-. 

V r max 

The average value of the transverse flow velocity, is very close to 

(2/3) (3f* x . 

The results are presented in Fig. 6, where we show the optimum values 
for the invariant time, r, and the flow parameter, /? max , together with their 
errors. If the results were consistent, then the fitted values would overlap. 
This is not the case. The fit to NA44 is far away from the fit to NA49, and 
the fit to STAR is far away from the fit to PHENIX. Amusingly, PHENIX is 
close to NA44, and STAR is close to NA49. The grand average of /? max from 
combined PHENIX and STAR is close to the average from combined NA44 
and NA49 experiments. These results indicate that even within a model 
capable of explaining the spectra one cannot conclude of a larger transverse 
flow at RHIC compared to SPS. 



5. Conclusions 

Here are our main points: 

1. The similarity of p^-spectra at SPS and RHIC in the data range sug- 
gests similar soft physics. We have argued that the combined present 
data do not lead to the conclusion of much larger transverse flow at 

4 We note that T and fis are the only independent thermodynamic parameters, since 
Us and fj,i follow from the conservation of strangeness and electric charge. 
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RHIC. Only the STAR data for p and K~ support this view, although 
with large experimental uncertainty. The PHENIX data are compat- 
ible with the same transverse flow as in SPS. 

2. The property of the spectra that there are more protons than 7r + at 
large p± can be seen not only in the RHIC data, but also in the SPS 
data. 



3. At the moment the experimental discrepancies between NA44 and 
NA49, and between STAR and PHENIX are of the same magnitude 
as the discrepancies between SPS and RHIC. With better data more 
accurate conclusion would be achieved. In particular, the RHIC mea- 
surements at lower energies will be very useful in verifying theoretical 
hypotheses. The use of models failing to reproduce the RHIC data 
should be, if possible, avoided in the modeling of detectors and in the 
analysis of the data. Also, the frequently made corrections for weak 
decays are not very useful, since these can be accounted for without 
difficulty in theoretical models. When comparing the corrected data 
much care is needed as to how the feeding from weak decays has been 
subtracted. 



4. We recall that the freeze-out temperatures used in the framework of 
thermal models have saturated (compare Eqs. (|l|) and @), which was 
not yet the case at AGS energies. Note, however, that the thermal 



approach works well for the AGS data [23, 24] and for elementary 



collisions f2q| . 



5. The transverse HBT radii and the slopes of the spectra are similar at 
SPS and RHIC, which implies similar transverse flow, or, in general, 
similar soft transverse physics. The particle yields increased by 50% 
from SPS to RHIC naturally allow for a 20% increase of the transverse 
size and flow. We also recall that the magnitude of the elliptic flow 
coefficient, V2, is similar from SPS to RHIC 2£j. 



6. The similarity of the soft physics at RHIC and SPS may be explained 
by the parton saturation phenomenon p7l |28|] (scenario (2) of Ref. 



[29]). If the onset of saturation occurs already at SPS energies, then 
the initial conditions for the multiparticle production are similar at 
SPS and RHIC, explaining the similarities discussed in this paper. 
The situation is reminiscent of the Hagedorn saturation (3^, [Jl], ^2| in 
elementary collisions, where the further increase of the collision energy 
does not lead to increased temperature. 
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The property of the saturation of soft physics, or at least the very weak 
dependence on the collision energy, should help to verify and constrain 
various models. The incomplete list of the most popular approaches 
and ideas includes: thermal models 



38], hydrodynamic models |39, 4C, 41 
transport theories El, and saturation models 
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Fig. 1. Comparison of the mid-rapidity pj_-spectra of ir , K , p, and A for the 
most central collisions of Pb + Pb at ^/snn = 17 GeV (NA44, NA49) and Au + Au 
at = 130 GeV (PHENIX, STAR). The STAR data for tt" and K~, and the 

NA49 data are preliminary. 
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Fig. 2. Same as Fig. 1 for ir + , K + , p, and A. 
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Fig. 3. Same as Fig. 1 for the rescaled spectra. 
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Fig. 4. Same as Fig. 2 for the rescaled spectra 
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Fig. 5. The contours of A\ 2 = 1 for various fits of the rescaled mid-rapidity 
p_L-spectra to the exponential form Aexp(— mj_/A). Optimum values are denoted 
by crosses. The numbers in parentheses are the values of the x 2 P er degree of 
freedom at the optimum. The data for NA49 and the STAR data for 7r~ and K~ 
are preliminary. 
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Fig. 6. The contours of A% 2 = 1 for fits of the combined mid-rapidity £>j_-spectra 
to the single- freeze- out model of Refs. Jl5|, ^2|. The parameters r and /3™ ax are 
defined in the text. The optimum values of the parameters are denoted by crosses. 
The numbers in parentheses are the values of the \ 2 P er degree of freedom at the 
optimum. The NA49 data and the STAR data for n~ and K~ are preliminary. 



